arXiv:1501.01837vl [physics.flu-dyn] 8 Jan 2015 


The effect of velocity boundary conditions on the heat transfer and flow topology in 

two-dimensional Rayleigh-Benard convection 

Erwin P. van der Poel 1 , Rodolfo Ostilla Monico 1 , Roberto Verzicco 1,2 and Detlef Lohse 1 
1 Department of Physics, Mesa-h Institute, and J. M. Burgers Centre for Fluid Dynamics, 

University of Twente, 7500 AE Enschede, The Netherlands 
2 Dipartimento di Inqeqneria Industriale, University of Rome u Tor Verqata”, Via del Politecnico 1, Roma 00133, Italy 

(Dated: January 9, 2015) 

The effect of various velocity boundary condition is studied in two-dimensional Rayleigh-Benard 
convection. Combinations of no-slip, stress-free and periodic boundary conditions are used on both 
the sidewalls and the horizontal plates. For the studied Rayleigh numbers Ra between 10 8 and 10 11 
the heat transport is lower for F = 0.33 than for F = 1 in case of no-slip sidewalls. This is surprisingly 
opposite for stress-free sidewalls, where the heat transport increases for lower aspect-ratio. In wider 
cells the aspect-ratio dependence is observed to disappear for Ra > 10 10 . Two distinct flow types 
with very different dynamics can be seen, mostly dependent on the plate velocity boundary condition, 
namely roll-like flow and horizontal zonal flow, which have a substantial effect on the dynamics and 
heat transport in the system. The predominantly horizontal zonal flow suppresses heat flux and is 
observed for stress-free and asymmetric plates. Low aspect-ratio periodic sidewall simulations with 
a no-slip boundary condition on the plates also exhibit zonal flow. In all the other cases, the flow 
is roll-like. In two-dimensional Rayleigh-Benard convection, the velocity boundary conditions thus 
have large implications on both roll-like and zonal flow that have to be taken into consideration 
before the boundary conditions are imposed. 

PACS numbers: 47.27.-i, 47.27.te 


I. INTRODUCTION 


Natural convection is commonly modelled by the 
Rayleigh-Benard (RB) system [Tf0] to reduce the com¬ 
plexity of the problem while preserving the most inter¬ 
esting features. This is a system in which a fluid is 
heated from below and cooled from above; gravity is act¬ 
ing in the plane normal to the energy-injecting or energy- 
extract inglower and upper plates. The temperature in¬ 
duced density differences result in a (temperature depen¬ 
dent) vertical forcing. Under the Boussinesq approxima¬ 
tion, in which the fluid properties are assumed to be in¬ 
dependent of temperature except for the buoyancy term, 
the governing equations in dimensionless form are: 


d t Ui + UjdjUi = -dip + \j ^d]ui + 0e 3 , 
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where £3 is the unit vector anti-parallel to gravity, iq is 
the velocity normalized by the freefall velocity ^/gflAL, 
t is the time normalized by the freefall time yjL/(g/3 A) 
and 0 is the temperature normalized by A, the tempera¬ 
ture difference between top and bottom plate, and shifted 
such that 0 is in the range 0 < 0 < 1. The control param¬ 
eters of the system are the non-dimensional temperature 
difference, i.e. the Rayleigh number Ra = gf3 AL 3 / {y k) , 
the Prandtl number of the fluid Pr = v/k and the aspect- 
ratio r = D/L, where L is the height of the sample and 
D its width, g the gravitational acceleration, /3 is the 
thermal expansion coefficient and v and k the kinematic 
viscosity and the thermal diffusivity, respectively. 


RB convection is used as a simple model for all types of 
natural convection, regardless of geometry. For 2D RB, 
at a moderate Ra = 10 s , the heat flux is expected to be 
independent of the aspect-ratio, and thus of the sidewall 
boundary conditions, beyond F >26 [5]. This value of 
T is expected to be lower for 3D RB [6|. Low aspect- 
ratio applications with a no-slip BC on the lateral walls, 
such as cooling fins and domestic heating, have moder¬ 
ate Rayleigh numbers (Ra < 10 10 ) that are accessible 
with contemporary direct numerical simulations (DNS). 
However, Nature is dominated by flows with large aspect- 
ratios and inaccessible high Ra [2 . 

In both experiments and numerics, the aspect-ratio 
T is commonly reduced to decrease the volume of the 
setup, as Ra depends only on the height of the system. 
In simulations, reducing the aspect-ratio for fixed Ra de¬ 
creases the computational complexity by (D(F 2 ) in three- 
dimensions, while in experiments it means a smaller setup 
volume. In both cases it makes higher Rayleigh numbers 
more accesible. Despite the additional effort required 
for large aspect-ratio, research in this direction has been 
done in simulations [6j and experiments CHS]. However, 
these endeavours were all limited in Ra due the large 
aspect-ratios used. In simulations one can approximate 
infinitely large F by using periodic boundary conditions 
at the sidewalls. This leads to the question how large 
the simulated domain has to be to sufficiently approxi¬ 
mate infinite T, as a truly infinite domain can in principle 
support any wavenumber while a periodic domain allows 
only for multiples of the fundamental one imposed by its 
size. This is addressed in this study in 2D RB, which is 
expected to require a larger domain than 3D and there- 
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FIG. 1. (Color online) Temperature field snapshots for different simulations. Red and blue indicate hot and cold fluid 
respectively. The colour varies between 0 — 0.2 and 0 — 0.8. Ra = 10 11 and Y — 1/2 for lateral periodicity, showing 
zonal flow (left panel), Ra = 10 10 and T = 0.33 for no-slip sidewalls, showing roll structures (center panel), and Ra = 10 10 and 
T = 0.33 for stress-free sidewalls (right panel). The plates are no-slip in all cases. Movies can be found in the supplementary 
material. 


fore is more demanding in terms of computational re¬ 
sources. 

Different temperature boundary conditions at the 
plates were studied by pm, where no significant difference 
between constant flux and isothermal plates was found in 
the higher Ra regime. Therefore, in this manuscript we 
focus only on the velocity boundary conditions. High 
Rayleigh number applications of RB usually do not have 
a no-slip lateral wall BC. Such lateral walls are unfortu¬ 
nately unavoidable in experiments due to flow confine¬ 
ment, but it is not required in numerics, although nu¬ 
merics tend to use no-slip sidewalls when replicating or 
complimenting experiments. Removing the no-slip BC 
on the lateral wall is actually beneficial in simulations, 
as it reduces computational requirements owing to the 
absence of the viscous BL and therefore no need to clus¬ 
ter grid points in the wall-normal direction. In any case, 
boundary conditions must be chosen in numerical stud¬ 
ies and one has to chose the boundary best suited for 
the intended study, but also the computational demand 
and possible introduced artifacts have to be taken into 
account. 

Despite its perceived importance, the effect of the ve¬ 
locity boundary conditions on the lateral walls has not 
been systematically studied for RB even though it is ex¬ 


pected to have a significant effect on the roll state, abso¬ 
lute heat flux and kinetic energy of the flow. In numerics, 
any expression which can be mathematically formulated 
can be used as a boundary condition, but here we focus 
on the three most physical options: no-slip (NS), stress- 
free (SF) and periodic (PD). For no-slip and stress-free 
sidewalls, an additional temperature BC has to be im¬ 
posed which in our case is adiabatic. 

Many recent studies on RB convection focus on the 
so-called “ultimate regime” QTHTjg, a regime of RB flow 
where the boundary layers are expected to become fully 
turbulent. It is expected that the scaling of integral quan¬ 
tities in this regime can be extrapolated to the very high 
Ra numbers present in geo- and astrophysical natural 
convection. In this manuscript, we address the question 
of how the lateral wall boundary conditions affect the flow 
and how we can extrapolate the results obtained to flows 
at arbitrary large driving, such as those in the ultimate 
regime. We find that, in essence, when one strives to 
model the high Ra applications of natural convection, the 
no-slip lateral wall BC and the low aspect-ratios of exper¬ 
iments and simulations may obscure phenomena which 
are seen in large aspect-ratio convection found in Na¬ 
ture. This research focuses on two-dimensional RB and 
it must be noted that three-dimensional RB might be- 
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FIG. 2. (Color online) Compensated Nu/Ra 2 ^ 7 as a function of Ra for stress-free and no-slip sidewalls and no-slip plates. The 
data is for aspect-ratios T =» 1 and T = 0.33. As expected, the stress-free Nu is consistently higher than for no-slip for the same 
value of T. For stress-free boundary conditions, the heat transport is higher for V = 0.33 than for T = 1 for all of the evaluated 
Ra, while for no-slip this is only observed at high Ra. 


have differently. In specific cases, qualitative arguments 
result in the expectation that three-dimensional RB does 
not show phenomena at all that are very strong in two- 
dimensional RB [16 j. 

One of these phenomena is zonal flow, which is a flow 
state with most of the kinetic energy stored in the hor¬ 
izontal motion pH CEE]. The topology consists of two 
shearing layers of fluid, one at the top and one at the bot¬ 
tom, that move in opposite horizontal directions. Zonal 
flow strongly suppresses heat flux, and correspondingly 
the heat flux as a function of time behaves very differ¬ 
ently. It can either be bursting or be sustained. In the 
first case quiescent periods of negligible convective heat 
flux are interrupted by large bursts, while in the sec¬ 
ond case the heat flux is more sustained. Even in the 
sustained case, the heat flux is much lower than would 
be expected in the more conventional roll-like flow state 
with identical control parameters. Zonal flow is observed 
with stress-free, asymmetric and no-slip plates, though 
for no-slip plates, the size of the period in laterally peri¬ 
odic simulations has to be small. 

The organization of the manuscript is as follows: First, 
the boundary conditions are formulated and their conse¬ 
quences on the numerical algorithm are evaluated. Then, 
simulations dominated by roll-like dynamics are pre¬ 
sented in both thin (T < 1) and wide (T > 1) cells. 
Here, we compare the heat flux in no-slip, stress-free and 
periodic sidewalls as function of Ra. We quantify the 
behaviour of the heat flux as function of T for periodic 
boundary conditions in order to find the aspect-ratio for 
which periodicity does not introduce any artifacts. Fi¬ 
nally, we study the effect of stress-free plates, both in 
an asymmetric setup with only one stress-free plate, and 
also with both plates having stress-free velocity bound¬ 
ary conditions. Here the heat flux as a function of time is 


dominated by large bursts with quiescent intervals. For 
all simulations a Prandtl number equal to unity is used, 
unless stated otherwise, to prevent the nesting of the ki¬ 
netic into the thermal boundary layer and vice versa, and 
to reduce the parameter space to a manageable size. 

II. NUMERICAL METHOD 

The direct numerical simulations (DNS) were per¬ 
formed with a two-dimensional version of a second or¬ 
der finite difference scheme |19| [20]. The underlying nu¬ 
merical scheme has been used for simulations of three- 
dimensional RB mm, Taylor-Couette 22 and other 
flows [23 . The two-dimensional version of the code was 
validated by comparison to a fourth order finite differ¬ 
ence code j24: that on its turn was compared with quasi 
two-dimensional experiments. Specifically, the Nu(Ra) 
data for T = 0.33 with no-slip boundary conditions can 
be compared one-to-one with the fourth order code and 
turn out to be equal within error [5 . In addition, results 
for no-slip boundary conditions in a T = 1 cell at Pr = 1 
agree well with the fourth order compact scheme of [25] . 

The Nusselt number Nu is calculated using 

Nu = v / RaPr {u z 0) Xtt - {d z 6) Xjt , (3) 

where x and z are the horizontal and vertical coordi¬ 
nate, respectively. Nu is a time average and one could 
denote it as (Nu)* and remove the time average from 
its definition, however Nu is used conform to RB liter¬ 
ature. For the instanteneous Nu, the notation Nu(£) is 
used. Nu is calculated as both an area average of the 
convective term Nu = V^RaPr (u z ) a, t + 1 and a line av¬ 
erage on both the horizontal plates of the diffusive term 
Nu = — (d z 0\ z =o : L)x,t- Note that the area average in 
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FIG. 3. (Color online) a) Nu as a function of T for stress-free sidewalls and no-slip plates. The data is for Ra = 10 8 . At 
r « 1 there is a transition region between a roll-state with one roll and a roll state with two horizontally aligned rolls, which is 
reflected by hysteretic jumps in Nu. b) Compensated (Nu — 1)/Ra 2 ^ 7 as a function of Ra for stress-free and periodic sidewalls 
and no-slip plates. Data is depicted for T = 1.5 and T = 2. The effect of the aspect-ratio on Nu disappears for high Ra while 
the effect of the sidewall boundary condition remains. 


two dimensions corresponds to a volume average in three 
dimensions. These Nu’s are verified by calculating Nu 
through the kinetic and thermal dissipation rates 0126]. 
A discrepancy between any of the altogether five differ¬ 
ent ways of calculating Nu indicates either insufficient 
statistical convergence or an underresolved simulation 
and so comparing these is a useful test. In addition, 
resolution checks, in which simulations are redundantly 
run using an increased resolution, are used to prove suffi¬ 
cient resolution for some cases. Additional details of the 
simulations can be found in the Appendix. 

The boundary conditions chosen are either no-slip, i.e. 
Ui = 0 at the walls, stress-free walls, i.e. iqn* = 0 , where 
rii is the vector normal to the wall, and d T Ui = 0 , where 
r is the direction tangent to the wall, and lateral period¬ 
icity, i.e. Ui(x = 0 , 2 :) = Ui(x = D,z). The temperature 
boundary condition on the plates are isothermal in all 
cases. For the no-slip and stress-free sidewalls, as ad¬ 
ditional temperature boundary condition we chose adia- 
baticity, i.e. d n 0 = 0 . 

The use of different boundary conditions has a signif¬ 
icant bearing on the algorithm to be chosen to solve the 
Boussinesq equations. The main difference is that for 
stress-free and periodic boundary conditions grids with 
uniformly spaced points can be used without loss of effi¬ 
ciency. In contrast, for no-slip, resolving the sharp gra¬ 
dients in the resulting boundary layers with a uniform 
grid-distribution will over-resolve the bulk. It is there¬ 
fore more efficient to use a non-uniform grid to cluster 
points near the boundary layers. However, solving an 
equation on a non-uniform grid is generally more costly 
in terms of computational requirements when compared 
to solving the same equation on an equally sized uniform 


grid. This is by virtue of the uniform grid that allows 
for the use of a Fourier transform to solve the pressure 
equation to impose the incompressibility of the flow. 


III. ROLL DOMINATED CONVECTION 

Thin cells 

For stress-free sidewalls, Nu for Y = 0.33 is higher 
than for T = 1 , which is in contrast with the no-slip case, 
where Nu for Y = 1 is higher than for Y = 0.33. It must 
be emphasized that this strong aspect-ratio dependence 
is a 2D effect not observed in 3D [ 6 ]. Even for V < 1 it is 
observed in 3D experiments that there is no difference in 
Nu between T = 1 and T = 1/2 cells for identical Ra and 
Pr [27l[28] , This is despite the fact that both in 2D and 
3D the amount of vertically stacked rolls changes with 
decreasing aspect-ratio M- 

The observation that the dependence of Nu on T dif¬ 
fers between no-slip and free-slip sidewall boundary con¬ 
ditions can be rationalized as follows. In the stress-free 
case, the aspect-ratio T = D/H can be interpreted as the 
ratio between the proportion of the enclosing wall that 
draws energy out of the system (no-slip plates) and the 
proportion that has no effect on the flow (stress-free side- 
wall). With decreasing T, this ratio decreases and thus 
there is less drag experienced by the flow at lower T, re¬ 
sulting in a higher Nu for stress-free sidewalls. In case of 
no-slip sidewalls, following the same reasoning, one ex¬ 
pects that the Nu does not depend on T provided that 
the flow structure remains the same. The difference in Nu 
for no-slip sidewall between thin and square (T = 1 ) cells 
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FIG. 4. (Color online) a) Nu vs T for periodic sidewalls and no-slip plates for Ra G {10 s , 10 9 ,10 10 ,10 11 }, see legend in b). 
Jumps in Nu can be seen around V « 1.25, T « 1.25, T « 1.10 and T w 1.05 for Ra = 10 s , Ra = 10 9 , Ra = 10 10 and Ra = 10 11 , 
respectively. The Nusselt number statistics at lower T than the jump are more difficult to converge than the other data points 
due to the bursting nature of Nu(t) and have a larger error. They are included to indicate the change of flow state. In b) Nu 
is compensated with Nu(T = 3). 


must then be a result of a changing flow state. While in 
T = 1 cells a single roll is present, in T = 0.33 cells there 
are four vertically stacked rolls, see figure [TJ d. In 2D, the 
heat flux through the system depends on the amount of 
vertically stacked rolls m , which explain the difference 
in Nu. 

There are exceptions to this trend, as there are specific 
cases in which an aspect-ratio reduction resulted in an 
increased heat flux for no-slip walls [30] . In addition, for 
the no-slip sidewall data seen in figure [2] Nu for T = 0.33 
is lower than Nu for T = 1 for most of the Ra range 
explored. However, at Ra = 4.64 x 10 10 this is opposite, 
which has been attributed to the break up of the large 
scale rolls in the flow [5]. 

We can further quantify the trend of Nu(T) for thin 
cells and stress-free sidewalls, and compare against the 
no-slip case. For no-slip sidewalls, Nu(T) shows relatively 
large jumps between vertically stacked roll states with Nu 
generally invariant with T (for T < 1) provided that the 
roll state does not change m • In figure [3^i the corre¬ 
sponding results for stress-free sidewalls can be seen. In 
contrast to the no-slip case, no jumps can be seen and 
Nu decreases monotonically with T up to T = 0.95. At 
T « 1 the flow transitions from a roll-state with one to a 
roll-state with two horizontally aligned rolls. There ap¬ 
pears to be multistability in this transition region. The 
absence of jumps for T < 1 can be directly attributed to 
the fact that for T <0.95 the system is exclusively in a 
single roll state while for T = 1 the system is in a roll 
state with two horizontally arranged rolls, resulting in a 
jump in Nu. 


Wide cells 

Small T cells are dominated by the presence of the side- 
wall. However, even if one could expect this dependence 
to progressively vanish for the case of large T, Nu for 2D 
RB systems with no-slip sidewalls has a strong T depen¬ 
dence up to very large values of T [5] (i.e. T < 26 for 
Ra = 10 8 ). Here we repeat the analysis of ref. [5], but 
now for cells with lateral periodicity, to check whether 
the T dependence is still present at such high values of T 
and study the scaling of Nu(Ra). 

In figure [ 3)3 Nu against Ra for laterally periodic and 
stress-free sidewalls are compared for T = 1.5 and T = 2. 
Nu for stress-free sidewalls is higher than for periodic 
sidewalls. This difference can be attributed to two ef¬ 
fects, namely the roll state imposed by the impermeable 
walls and the drag from the sidewall. The amount of 
horizontally stacked rolls for periodic sidewalls, in con¬ 
trast to impenetrable sidewalls (NS and SF sidewalls), 
is restricted to even numbers. At T = 2, for the stud¬ 
ied parameter range, the periodic simulations showed two 
horizontally stacked rolls while the stress-free simulations 
showed three rolls. For both simulations, the initial con¬ 
ditions consisted of two rolls, to ensure that the heat flux 
is not affected by a different roll state. 

Nu heavily depends on the amount of horizontally 
stacked rolls m and this effect most likely contributed to 
the difference in Nu between the evaluated cases. How¬ 
ever, there could also be other causes of the difference 
in Nu. In order to find out, T = 1.5 was simulated for 
the same Ra as for T = 2. At this aspect-ratio both lat¬ 
erally periodic and stress-free sidewall simulations have 
a roll state consisting of two rolls. It can be seen that 
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FIG. 5. (Color online) Temperature field snapshots for Ra = 2.15 x 10 8 and r = 1.5 with no-slip plates and stress-free sidewall 
boundary conditions of two different statistically steady states. The colormap ranges from 0 < 6 < 1. The area averaged mean 
temperature (0)A,t is 0.44 (left panel) and 0.56 (right panel). 


the difference in Nu between PD and SF is smaller, as 
expected, and that the difference in stress at the bound¬ 
aries alone results in a lower Nu for laterally periodic cells 
as compared to stress-free sidewalls. 

The figure shows that Nu becomes T-independent for 
both BC’s at sufficiently large Ra. This finding adds the 
conclusion that not only laterally periodic cells lose the 
T dependence, but that this is also valid for stress-free 
sidewalls at T « 2 for Ra > 10 10 . We now focus on the 
aspect-ratio dependence of periodic sidewall simulations 
at different Ra. 

Periodic boundary conditions are used to approximate 
infinite aspect-ratio, i.e. an absence of finite size effects, 
on a finite size domain. Thus, the aspect-ratio in periodic 
sidewall simulations does not reflect a physical domain 
but instead the period size of the computational domain. 
This aspect-ratio cannot be arbitrarily small due to two 
reasons: i) there must always be an even number of hor¬ 
izontal rolls in a periodic system, and ii) wavelengths 
larger than the domain are not simulated and thus not 
accounted for. However, as the aspect-ratio is propor¬ 
tional to the computational cost, it is of interest to study 
the minimum aspect-ratio r m i n that still produces phys¬ 
ical solutions. When considering three dimensional (3D) 
flows with a large scale circulation (LSC), periodicity can 
be imposed in two distinct dimensions. The periodic di¬ 
mensions can either be parallel or perpendicular to the 
vorticity vector of the LSC. It is expected that parallel 
periodicity has less effect on the flow than perpendicular 
periodicity. Parallel periodicity is often used in Taylor- 
Couette simulations [22] . where, unlike in 3D RB [3T], the 
orientation of the LSC is known a priori. In 2D RB this 
periodicity is unavoidably imposed in the perpendicular 
direction and thus a strong effect is expected. 

In figure [4] Nu as function of the period T is shown 


for four different Ra. One would expect that r m i n « 2 
as the minimum number of rolls is two and the aspect- 
ratio of one roll is approximately unity [5]. It is therefore 
surprising that a weak dependence of Nu, i.e. < 5%, on 
T is seen from r min significantly lower than two, where 
T m in is indicated by the jump in Nu(Ra). For Ra = 10 9 
the aspect-ratio dependence is larger than for higher Ra. 
This is in particular visible in figure IIb At Ra = 10 9 , 
Nu differs approximately 8% between T = 1.5 and T = 2, 
while this difference vanishes at Ra = 10 10 . r m i n appears 
to decrease for increasing Ra. We expect this to be a 
result of a weaker LSC m and decreased structure size 
at higher Ra. 

We wish to highlight a peculiarity observed for stress- 
free sidewall and no-slip plates. Namely, that the mean 
temperature can be (0)A,t 7 ^ 0.5. This can be seen in 
more detail in figure [5] Even though there is top bot¬ 
tom symmetry in the boundary conditions and in the 
equations there is a deviation from the expected value 
(0)A,t = 0.5. As the system is not turbulent enough, 
it might take infinite time for the system to explore all 
possible states and the symmetry can be sponteneously 
broken due to the different flow states. The system is 
in a horizontally stacked double roll state where the hot 
plumes are moving up at the center of the no-slip plate 
and the cold plumes move down along the stress-free side- 
walls. In this case ( 0)A,t > 0.5, as the hot plumes are 
well mixed by turbulence while the cold plumes seem to 
lose their heat through molecular diffusion, see figure [5^i. 
A similar but opposite effect can be obtained in other 
realizations of the same system with different initial con¬ 
ditions. In steady state axisymmetric convection with 
stress-free plates and sidewalls [32] the core temperature 
is smaller than 0.5 in case of a rising central plume. 

We have checked the dependence of the deviation of the 













7 




x/L 


x/L 


FIG. 6. (Color online) Temperature field snapshots with velocity vectors superimposed for stress-free plates and periodic 
sidewalls for Pr = 1 and at Ra = 10 6 (left panel) and Ra = 10 8 (right panel) in a T — 2 cell. The temperature ranges 
from 0 < 0 < 1. The size of the arrows is proportional to the absolute velocity. The flow is roll-like in a) and zonal in b). 
Corresponding movies can be found in the supplemental material. 


mean temperature from 0.5 on the central plume location 
by performing the transformation Ui(x, z) uflx , L — z) 
and 0(x,z) 1 — 0(x,L — z) on the hot solution of 

T = 1.5, Ra = 2.15 x 10 8 . This transformation flips 
the velocity and temperature fields upside-down and in¬ 
verts the temperature. While before the transformation 
in statistically steady state (0)A,t = 0.54 > 0.5, converg¬ 
ing the simulation with the transformed field results in 
(0)A,t — 0.46 < 0.5, of which a snapshot can be found in 
figure [ 5 J 3 . This indicates that the deviation from 0.5 is 
a flow state effect heavily dependent on the initial con¬ 
dition. Even though the time and volume averaged heat 
flux through the system is identical in both cases, the 
area-averaged and local temperature strongly differ be¬ 
tween these two flow states. 


IV. ZONAL FLOW 
Stress-free plates 

Up to now we have focused on the case with no-slip 
plates and varied the sidewall boundary conditions. Of 
course, another degree of freedom of the boundary condi¬ 
tions are the horizontal plates. Using a stress-free bound¬ 
ary condition instead of no-slip results in huge differences 
in both flow topology and heat flux. Stress-free plates are 
used in theory to simplify boundary conditions, as well 
as infinite (periodic) systems. This BC is a critical re¬ 
quirement for theorems that use the vanishing vorticity 
condition at the plates, which goes all the way back to 
original analysis by Rayleigh, but it is also used for more 
recent work, such as the finite-Pr upper bound theorem 


of [33] that gives Nu < 0.2891Ra 5 ^ 12 . 

We performed simulations with (horizontal) stress-free 
plates and used lateral periodicity with T = 2 to ap¬ 
proximate an infinite system. For Ra > 10 7 the flow 
topology consists of two shearing layers with predomi¬ 
nantly horizontal motion. The flow at the lower half of 
the domain moves in opposite direction as the top. This 
is because the total momentum must be conserved due to 
the absence of momentum transport to the boundaries. 
The flow is in a state for which Nu shows bursts with 
long “quiet” intervals where Nu(£) ~ 1 . This flow state 
is referred to as zonal flow pa m mi mg, similar to 
the zero-wavenumber flows observed in toroidal plasmas 
and planetary atmospheres, e.g. on Jupiter. The strong 
horizontal motion is driven by the so-called “shearing 
mechanism”. Plumes emitted from the one plate move 
vertically at first until they are deviated horizontally by 
the opposite plate, driving the horizontal motion. 

Depending on Pr, there is a threshold Ra beyond which 
the flow transitions from roll-like to zonal flow. Hystere¬ 
sis is observed at this transition [18]. In figure [ 6 ] snap¬ 
shots can be seen of the flow below and above this thresh¬ 
old at Pr = 1 . At Ra = 10 6 the flow is still roll-like; a 
large scale circulation with accompanying plumes. In¬ 
creasing Ra to Ra = 10 8 the flow transitions to a zonal 
flow. In the corresponding snapshot in figure [ 6)3 the flow 
is in the Nu(£) ~ 1 interval between bursts, in which a 
linear temperature gradient can be seen with strong hor¬ 
izontal flow. This interval phase is interrupted by large 
bursts in Nu(t). During these burst phases there is high 
mixing before the interval phase is re-initiated. Zonal 
flow is only observed in two dimensions. In our three- 
dimensional simulations using SF plates, periodic side- 
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FIG. 7. (Color online) a) Nu as a function of dimensionless time in freefall time units at T = 1/2, for Ra = 1 X 10 8 and 
Ra = 1 x 10 11 , NS plates and PD sidewalls. At Ra = 10 8 we observe bursting behaviour while at Ra = 10 11 we do not: The 
zonal flow is bursting for lower Ra and sustained for higher Ra. Note the large difference between the time scale of the bursts 
at Ra = 10 8 and the fluctuations at Ra = 10 11 . b) Probability density function of the time interval between the bursts At for 
Ra — 10 7 and Ra — 10 8 . 


wall and a cartesian geometry with T = 1, no zonal-flow 
has been observed for Ra G {10 7 ,10 8 ,10 9 }. Bursts in 
Nu have been observed for wind-reversals [36] although 
in that case the bursts are not preceded by “quiet” i.e. 
Nu ^ 1 flow states. 

The zonal flow state is detrimental for obtaining 
Nu(Ra) scaling as Nu is close to one and, due to the 
intermittent nature of the bursts, time-averaged quanti¬ 
ties are very hard to converge, requiring simulations of 
several tens of thousands of large eddy turnover times. In 
order to obtain a converged Nu to see how the system ap¬ 
proaches the theoretical upper bounds, we prevent zonal 
flow by increasing Pr as this increases the shear between 
the bottom and top layer. Increasing to Pr = 40 pre¬ 
vents zonal flow for Ra < 10 10 . In the low Ra regime 
the flow is composed of two rolls with very thin plumes 
streaks in between, very similar to infinite Pr flows with 
SF plates [37] . The absolute Nu for a given Ra is nearly 
double of that in the case of no-slip plates. However, the 
Nu(Ra) scaling exponent is similar, namely 0.28, much 
smaller than the upper bound exponent of 5/12 « 0.42 
[3]J], rigorously derived for these boundary conditions. 


Suppressed periodicity 

Zonal flow is not inherent to the stress-free plate 
boundary condition. It is also observed for no-slip plates 
when T is too low for the lateral periodicity and the flow 
becomes dominated by artifacts. For T < r m i n the flow 
is zonal for all the studied Ra. First, we focus on T = 1, 
which is close to r m ^ n where the flow transitions from 
zonal to roll-like, see figure [4] A strong characteristic of 
zonal-flow as compared to roll-like flow is that the mean 
horizontal velocity {u x ) Xy t as a function of z/L is non¬ 
zero m- This can be seen in figure [8 ^l for T = 1. The 


two shearing layers are reflected in (u x ) Xj t by the two 
peaks with opposite velocity. The velocity profile seems 
nearly independent of Ra. Unlike with stress-free plates 
m, the velocity in this case vanishes at the plates, in¬ 
herent to the no-slip boundary condition. It thus appears 
that zonal flow is not intrinsic to stress-free plates. When 
the lateral period is too small to allow roll-like flow, the 
sideways deflected plumes have no freedom to return to 
their original plates. The plumes therefore have to move 
sideways, driving a zonal flow. 

Decreasing the aspect-ratio to T = 0.5 increasingly 
suppresses the flow. In figure [7^, Nu(t) for two different 
Ra is depicted. Nu(t) for Ra = 10 8 shows the Nu bursts 
with a peak value of approximately Nu(t pea k ) ~ 50. For 
Ra = 10 11 the features of Nu(t) are completely different. 
Here, large symmetric fluctuations around the mean can 
be seen, characteristic for common RB systems. A sim¬ 
ilar transition is seen for zonal flow between stress-free 
plates m- However, there the zonal flow type transitions 
from bursting to sustained when Pr is increased. Here, 
the transistion between the bursting and the sustained 
Nu is at Ra = 10 10 for constant Pr = 1 and T = 0.5, 
with bursts found at lower Ra. At this same Ra = 10 10 
the vertical profile of (u x ) Xj t starts to differ from the pro¬ 
files at lower Ra (figure For lower Ra, the horizontal 
motion is only in one direction, with a maximum at mid¬ 
height. At higher Ra, the velocity near the plates opposes 
the motion at mid-height and three regions of different 
horizontal motion are present. This shows that also for 
zonal flow there are different states, strongly dependent 
on the boundary conditions and correspondingly T. 

The scale of the buoyant thermal plumes decreases for 
both increasing Ra and increasing Pr, allowing for more 
freedom compared to the fixed periodic length scale. In 
the T = 0.5 case, increasing Ra prevents the bursting 
flow state. In the sustained zonal flow state, there are no 
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FIG. 8. (Color online) Mean horizontal velocity (u x ) Xl t, normalized by its maximum, as a function of vertical position z/L for 
different Ra (see legend) in a) F = 1 cells with lateral periodicity and no-slip plates, b) T = 0.5 cells with lateral periodicity 
and no-slip plates and c) T = 2 cells with lateral periodicity and asymmetric plates. 


bursts in Nu. Most likely it is the heat buildup during 
the stagnant phase that results in the bursts. This is 
supported by the fact that for larger T = 1, the zonal 
flow is additionally sustained. The flowfield snapshot of 
Ra = 10 11 is displayed in figure [lji. The flow topology in 
this snapshot is very different from typical RB (see e.g. 
figure HR even though Nu(t) has a similar structure, 
its mean value is substantially lower, indicating that the 
suppression of heat flux still is significant at higher Ra. 

The interval between the bursts strongly depends on 
Ra. In figure [TJd, the probability density functions for 
Ra = 10 7 and Ra = 10 8 can be seen. The interval du¬ 
ration At is extracted using a peak finding algorithm 
that finds local maxima. The minimum peak height and 
minimum peak distance are tuned by hand for each Ra. 
For increasing Ra, the mean interval increases while the 
variation appears approximately constant. For Ra > 10 8 
the mean At becomes so large that obtaining statisti¬ 
cal convergence requires huge amounts of CPU time. It 
therefore remains unknown how the mean time interval 
between bursts increases up to the point that the burst¬ 
ing behaviour disappears. In addition no converged Nu 
can be obtained as a result of the increasing time interval. 

These low T PD sidewall simulations are the only cases 
we encountered with zonal flow when both plates are no¬ 
slip in two-dimensions. However, it has been observed in 
three-dimensional RB experiments in an annular cylin¬ 
der m, where all the enclosing boundary conditions are 
no-slip with periodicity introduced in the azimuthal di¬ 
rection. 


Asymmetric plate boundary conditions 

The final variation we consider on the plates is using 
asymmetric velocity boundary conditions. These con¬ 
ditions are rarely studied, but have many applications; 
the ocean or a cooking pan can be thought of as having 
stress-free top at their free surfaces and a no-slip bot¬ 


tom. Asymmetric boundary conditions break the top- 
down symmetry in RB even in the Boussinesq approxi¬ 
mation. In conjunction to asymmetric boundary condi¬ 
tions at the plates, we take a laterally periodic system, 
to avoid as much as possible sidewall effects. 

Snapshots of the flowfield are displayed in figure [9^i. 
Here, it can be seen that we obtain partly zonal flow 
near the top plate where the BC is stress-free. Indeed, 
the system displays similar behaviour to the case where 
both plates are stress-free, namely bursts in Nu with 
Nu « 1 intervals. The main difference between stress- 
free BCs and the asymmetric setup is obviously near the 
no-slip bottom plate. In symmetric zonal flow, the mech¬ 
anism that drives zonal flow is conceptually understood 
as plumes from both plates driving the zonal flow at the 
opposite plate. However, the zonal flow present in this 
asymmetric setup has no plumes originating from the 
bottom plate, see[9^. This shows that even without these 
plumes, the zonal flow in the upper, stress-free, layer can 
be driven by the interaction between the stationary layer 
at the bottom and the descending cold plumes alone. 

In our simulations, we see this mixed zonal-flow state 
for Pr = 1 and Ra > 10 9 . For lower Ra, the system 
displays large scale circulations without zonal flow, re¬ 
sulting in a more typical behaviour of Nu(t) fluctuat¬ 
ing around the mean without periods where Nu ~ 1. 
The thermal plumes in this Ra regime have a different 
shape; the plumes originating from the no-slip plate are 
surrounded with high vorticity while the others are not. 
The plumes that originated at the stress-free plate ap¬ 
pear like long streaks that are less well mixed and have 
sharper gradients compared to their NS counterparts. 

Figure [9 ]d shows Nu as a function of Ra. As these Nus- 
selt number data are very hard to converge the typical 
averaging time is at least one order of magnitude larger 
than for roll-like flows. The difference between Nu with 
and without zonal flow can be clearly seen. Moving up¬ 
wards in Ra, i.e. using the solution field from a lower Ra 
simulation as an initial condition, the flow transitions at 

























10 



FIG. 9. (Color online) a) Temperature field snapshot for asymmetric plates BCs and periodic sidewalls at Ra = 4.64 x 10 8 . The 
top boundary is stress-free and the bottom boundary is no-slip. . Movies of roll-like and zonal flow with asymmetric BCs can 
be found in the supplemental material, b) Nu as a function of Ra for asymmetric plate BC’s and periodic sidewalls. Simulations 
where the initial conditions is taken from the lower Ra neighbouring data point are indicated with rightward-pointing triangles 
and where the initial condition is taken from higher Ra are indicated with leftward-pointing triangles. 


Ra = 10 9 . Moving backwards, the transition back to¬ 
wards a flow without zonal flow is at Ra = 10 7 . This 
difference indicates that there is hysteretic behaviour at 
this transition, similar to the double stress-free plates 
case m- The vertical profile of (u x ) Xj t can be seen in 
figure for the cases with asymmetric plates that have 
zonal flow. Below the forward transition at Ra = 10 9 , 
the horizontal motion is exclusively in one direction. For 
higher Ra, similar to the case described in the previous 
chapter, a layer is observed near the no-slip plate that 
has an opposite horizontal motion. At higher Ra, the 
plumes coming from the stress-free cold plate reach the 
no-slip plate and drive the zonal flow in the bottom of 
the cell, while at lower Ra the plumes do not reach that 
far into the flow and are dissipated before they can drive 
the zonal flow. 

As the up-down symmetry is broken it is expected that 
the area and time averaged temperature {0)a, t can devi¬ 
ate from the value of 0.5 commonly found in symmetric 
RB [39], similar to what is observed in non-Oberbeck- 
Boussinesq RB convection [38] . This quantity is plotted 
as a function of Ra in figure[T0^, where it can be seen that 
also within the Oberbeck-Boussinesq approximation the 
deviation from 0.5 can be substantial, as seen here with 
asymetric plate boundary conditions. In addition, there 
is no clear dependence on Ra except for a large jump 
that corresponds to the transition between roll-like and 
zonal flow. It appears that for roll-like flow at lower Ra, 
the plate with the stress-free BC (in this case the cold 
plate) has a dominant effect on the bulk temperature, 
as (0)A,t < 0-5- The absence of the viscous boundary 
layer on this plate must therefore have a stronger effect 
on {0) a, t than the increased vorticity at the no-slip plate. 
At the stress-free plate, the full thermal boundary layer is 
exposed to the LSC while at the no-slip plate it is (partly) 
nested in the viscous boundary layer. In contrast, in case 


of zonal-flow at higher Ra, there is a negligible deviation 
of the mean temperature from the symmetrical value of 
0.5. 

In figure [Tob the vertical temperature profiles (0) Xjt as 
function oizjL can be seen for various Ra. The no-slip 
and stress-free boundary layers are located near z/L = 0 
and z/L = 1, respectively. For Ra = 10 7 and Ra = 
10 8 the bulk profile is clearly colder than 0.5. The two 
different boundary layer profiles are significantly different 
from each other for both these Ra. The no-slip plate 
boundary layer is less steep and larger than the stress- 
free one. A small overshoot in 9 can be seen for Ra = 10 7 , 
which is typical for low Ra RB and not too small Pr. The 
boundary layer profile for Ra = 10 9 deviates strongly 
from typical RB BL profiles m ■ It must be noted that 
this is a time average and so this profile is affected by 
both flow phases; the zonal flow phase and the burst. As 
there is negligible vertical motion in the zonal flow region, 
a near linear temperature profile can be seen, confirming 
that the convective heat flux is negligible and that the 
heat flux is dominated by molecular diffusion. 


V. SUMMARY AND OUTLOOK 

As was seen in the previous sections, the veloc¬ 
ity boundary conditions in two-dimensional Rayleigh- 
Benard convection have a large effect on the flow topol¬ 
ogy and consequently, the Nusselt number. This is in con¬ 
trast with the temperature boundary conditions, which 
have been shown to have neglible effects when switching 
between isothermal and constant flux boundary condi¬ 
tions on the plates m ■ 

The flow can be classified as roll-like and zonal. The 
roll-like flow state is composed of a large scale circulation 
interacting with thermal plumes. The resulting Nu(t) 
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FIG. 10. (Color online) a) The area -and time averaged temperature (0)A,t as a function of Ra for asymmetric plate boundary 
conditions. Only the runs where the initial condition is taken from a lower Ra are included, b) The horizontal -and time 
average temperature {0) x ,t as a function of the vertical coordinate zjL for Ra = 10 7 , Ra = 10 8 and Ra = 10 9 . The Ra = 10 9 
data in this plot is taken from a simulation with zonal flow while the others are in a roll-like state. The dashed black line in 
indicates {9) Xi t = 0.5. For both panels, the asymmetry is introduced by making bottom plate (i.e. zfL — 0) no-slip, while the 
top plate (i.e. zjL — 1) is stress-free. 


fluctuates around its mean. This flow state is most com¬ 
monly observed when the horizontal plates are subject to 
a no-slip boundary condition, which is the case in exper¬ 
iments. The zonal flow state is however very different. 
When looking at the Nu time series, we see large bursts 
of Nu surrounded by “quiet” intervals where Nu(t) « 1. 
During these intervals, vertical motion is absent and the 
flow largely moves horizontally. This is allowed by the use 
of lateral periodicity. Zonal flow is observed when either 
one or both plates have a stress-free boundary condition. 
Zonal flow is additionally observed for no-slip plates in 
a special case, i.e. when the sidewall is periodic and 
T < T min , and thus the flow has limited freedom. This 
results in a similar Nu(t) structure. A difference between 
these cases is that zonal flow is observed upwards from a 
critical Ra, while for the low T case with lateral periodic¬ 
ity, the bursting behaviour of Nu(t) is only observed for 
a Ra lower than a critical Ra. In all cases, zonal flow is 
only observed for periodic sidewalls in two-dimensional 
RB. 

We have quantified the effect of the sidewall on Nu 
for no-slip plates when the flow is roll-like using for no¬ 
slip sidewalls, stress-free sidewalls and lateral periodicity. 
As expected, Nu is higher for stress-free than for no-slip 
sidewalls over the evaluated Ra range. Surprisingly, in 
the stress-free case, Nu is higher for lower T, while the 
opposite happens for no-slip. This can be attributed to 
the absence of vertical roll states for stress-free sidewalls. 
In addition, stress-free boundary conditions result in heat 
flux, when compared to lateral periodicity for both T = 
1.5 and T = 2 and for both of these boundary conditions 
the aspect-ratio dependence of Nu disappears at high Ra. 

From analysing Nu(T) for lateral periodicity and no¬ 


slip plates, we can conclude that below a certain r m i n the 
flow becomes zonal. Roll-like flow is supressed, resulting 
in a lower Nu than larger T > r m i n . The type of zonal 
flow that results depends strongly on T, as zonal flows 
with two and three shearing layers are observed for T = 
1 and T = 0.5, respectively. r m i n decreases for higher 
Ra, most likely due to the decreased viscous stresses on 
the thermal plumes. This allows the plumes to require 
less freedom to move alongside each other in opposite 
directions, allowing for roll-like flow. 

Finally, in the case of asymmetric plate boundary con¬ 
ditions, where one horizontal plate has a no-slip bound¬ 
ary condition and the other plate a stress-free boundary 
condition, zonal flow can be observed from a critical Ra 
onwards, as was the case for two stress-free plates. Be¬ 
low this Ra, the flow remains roll-like, although the mean 
temperature is not 1/2 as the up-down symmetry is bro¬ 
ken due to the different velocity boundary conditions. 
The plate with the stress-free boundary condition has a 
larger effect on the mean temperature and it therefore de¬ 
viates from 1/2 towards the temperature of the SF plate. 
This same effect was also seen previously with stress- 
free sidewalls, even though the plate boundary conditions 
were equal. We wish to highlight that even though the 
effect is the same, they are due to very different mecha¬ 
nisms. The resulting vertical temperature profile is also 
completely different between the zonal flow and the roll¬ 
like flow regimes. In addition, we observed hysteretic be¬ 
haviour around the transition from roll-like to zonal flow, 
similar to the case with two stress-free plates, reported 

in [IE]. 

Choosing the boundary conditions for a simulation can 
be an involved, but very relevant, consideration, as re- 
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suits have shown here. Based on computational perfor¬ 
mance, the laterally periodic simulations with stress-free 
plates would be the first choice; removing most bound¬ 
ary layers. In two-dimensions this can result in zonal 
flow with very little resemblance to roll-like RB convec¬ 
tion, commonly observed in experiments. Zonal flow is 
avoided by introducing drag through the no-slip plate 
boundary condition. However, the roll-like convection 
that corresponds to this set of boundary conditions re¬ 
quires a minimum aspect-ratio, otherwise zonal flow will 
emerge. Above this, the aspect-ratio dependence, but 
not the boundary condition dependence, of the heat flux 
becomes negligible for Ra > 10 10 at least for stress-free 
and periodic sidewall. Introducing asymmetry in the sys¬ 
tem, between either the plates or the sidewall and the 
plates, can result in zonal flow but also surprising ef¬ 
fects as unexpected mean temperatures deviating from 
1/2. Most, if not all, of these boundary conditions ef¬ 
fects are qualitatively expected to be non-existent or at 
least much weaker in three-dimensions as compared to 
two-dimensions. 
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TABLE I. The columns from left to right indicate Ra, Pr, T, the wall BC, the plate BC, the resolution in horizontal and vertical 
direction N x x N z , the number of grid points in the thermal boundary layer #{n | n £ {1, 2, N z } A Z(n) < L/(2Nu)} 
(where Z(n) is the vertical grid point distribution), the average Kolmogorov length scale in the flow compared to the largest 
grid length used somewhere in the grid max(S x ,S z )/rj withr]/L = Pr 1 ^ 2 /[Ra(Nu — l)] 1 ^ 4 , the mean Nu over the full simulation 
and the averaging time in free fall time units t/. 


Ra 

Pr 

r 

wall BC 

plate BC 

N x x N z 

Nbl 

max(5 x ,5 Z ) 

Nu 

T f 

V 

10 Y 

1 

2.00 

PD 

AS 

512 x 256 

11 

0.47 

21.6 

750 

2.15 x 10 7 

1 

2.00 

PD 

AS 

512 x 256 

10 

0.59 

25.8 

900 

4.64 x 10 7 

1 

2.00 

PD 

AS 

512 x 256 

9 

0.76 

32.0 

900 

10 s 

1 

2.00 

PD 

AS 

512 x 256 

8 

0.95 

37.0 

900 

2.15 x 10 8 

1 

2.00 

PD 

AS 

1024 x 512 

15 

0.60 

43.0 

900 

4.64 x 10 8 

1 

2.00 

PD 

AS 

1024 x 512 

13 

0.76 

52.3 

1000 

10 9 

1 

2.00 

PD 

AS 

3072 x 1024 

118 

0.20 

9.8 

10000 

2.15 x 10 9 

1 

2.00 

PD 

AS 

3072 x 1024 

128 

0.25 

8.6 

14000 

10 9 

1 

1.00 

PD 

NS 

512 x 512 

19 

0.79 

27.9 

1700 

10 9 

1 

1.20 

PD 

NS 

512 x 512 

18 

0.96 

29.1 

1500 

10 9 

1 

1.25 

PD 

NS 

512 x 512 

16 

1.18 

54.8 

600 

10 9 

1 

1.30 

PD 

NS 

512 x 512 

17 

1.22 

53.7 

400 

10 9 

1 

1.40 

PD 

NS 

1024 x 512 

11 

0.66 

53.7 

1500 

10 9 

1 

1.60 

PD 

NS 

1024 x 512 

11 

0.74 

51.0 

1500 

10 9 

1 

1.80 

PD 

NS 

1024 x 512 

12 

0.82 

48.5 

1500 

10 9 

1 

2.00 

PD 

NS 

1024 x 512 

12 

0.91 

48.3 

1200 

10 9 

1 

2.20 

PD 

NS 

1024 x 512 

12 

1.00 

47.9 

1200 

10 9 

1 

2.40 

PD 

NS 

1600 x 512 

12 

0.70 

47.6 

900 

10 9 

1 

2.60 

PD 

NS 

1600 x 512 

11 

0.78 

54.8 

1000 

10 9 

1 

2.80 

PD 

NS 

1600 x 512 

11 

0.84 

53.4 

1300 

10 9 

1 

3.00 

PD 

NS 

1600 x 512 

11 

0.88 

50.9 

1200 

10 10 

1 

1.00 

PD 

NS 

1024 x 1024 

28 

0.85 

59.1 

400 

10 10 

1 

1.05 

PD 

NS 

1024 x 1024 

51 

0.86 

49.2 

500 

10 10 

1 

1.10 

PD 

NS 

1024 x 1024 

32 

1.08 

102 

1100 

10 10 

1 

1.20 

PD 

NS 

1024 x 1024 

18 

1.17 

99 

600 

10 10 

1 

1.40 

PD 

NS 

1440 x 1024 

18 

0.96 

95.8 

600 

10 10 

1 

1.60 

PD 

NS 

1440 x 1024 

18 

1.10 

95.7 

600 

10 10 

1 

1.80 

PD 

NS 

1440 x 1024 

18 

1.24 

96.4 

700 

10 10 

1 

2.00 

PD 

NS 

1440 x 1024 

42 

1.37 

94.9 

400 

10 10 

1 

2.20 

PD 

NS 

2048 x 1024 

19 

1.05 

93.7 

1000 

10 10 

1 

2.40 

PD 

NS 

2048 x 1024 

19 

1.15 

95.3 

1000 

10 10 

1 

2.60 

PD 

NS 

2048 x 1024 

19 

1.24 

93.1 

1100 

10 10 

1 

2.80 

PD 

NS 

2048 x 1024 

19 

1.34 

93.3 

1000 

10 10 

1 

3.00 

PD 

NS 

2048 x 1024 

23 

1.43 

92.6 

1000 

10 11 

1 

1.00 

PD 

NS 

2048 x 2048 

61 

0.82 

77.9 

300 

10 11 

1 

1.05 

PD 

NS 

2400 x 2048 

50 

0.92 

198 

200 

10 11 

1 

1.10 

PD 

NS 

2400 x 2048 

51 

0.95 

189 

200 

10 11 

1 

1.40 

PD 

NS 

2400 x 2048 

29 

1.23 

197 

300 

10 11 

1 

1.80 

PD 

NS 

2560 x 2048 

30 

1.48 

196 

200 

10 11 

1 

2.20 

PD 

NS 

4096 x 2048 

31 

1.11 

187 

200 

10 11 

1 

2.60 

PD 

NS 

4096 x 2048 

31 

1.32 

188 

200 

10 11 

1 

3.00 

PD 

NS 

4096 x 2048 

31 

1.52 

186 

200 

10 s 

1 

1.50 

PD 

NS 

512 x 256 

15 

0.66 

26.1 

600 

2.15 x 10 8 

1 

1.50 

PD 

NS 

1024 x 512 

24 

0.42 

33.1 

400 

4.46 x 10 8 

1 

1.50 

PD 

NS 

1024 x 512 

20 

0.54 

41.1 

400 

10 9 

1 

1.50 

PD 

NS 

1024 x 512 

17 

0.69 

51.3 

400 

2.15 x 10 9 

1 

1.50 

PD 

NS 

2048 x 1024 

31 

0.44 

63.3 

400 

4.64 x 10 9 

1 

1.50 

PD 

NS 

2048 x 1024 

27 

0.56 

75.5 

400 

10 10 

1 

1.50 

PD 

NS 

3072 x 1536 

38 

0.48 

94.4 

300 

10 s 

1 

2.00 

PD 

NS 

512 x 256 

15 

0.66 

26.1 

600 

2.15 x 10 8 

1 

2.00 

PD 

NS 

512 x 256 

9 

1.11 

31.2 

2600 

4.46 x 10 8 

1 

2.00 

PD 

NS 

512 x 256 

8 

1.42 

38.9 

2400 

2.15 x 10 9 

1 

2.00 

PD 

NS 

1024 x 512 

10 

1.17 

61.1 

1100 

4.64 x 10 9 

1 

2.00 

PD 

NS 

1024 x 512 

10 

1.50 

76.3 

1000 
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TABLE II. The columns from left to right indicate Ra, Pr, T, the wall BC, the plate BC, the resolution in horizontal and 
vertical direction N x x N z , the number of grid points in the thermal boundary layer #{n | n £ {1, 2, N z } A Z(n) < L/ (2Nu)} 
(where Z(n) is the vertical grid point distribution), the average Kolmogorov length scale in the flow compared to the largest 
grid length used somewhere in the grid max^, S z )/rj with r\/L — Pr 1//2 /[Ra(Nu — l)] 1 ^ 4 , the mean Nu over the full simulation 
and the averaging time in free fall time units t/. 


Ra 

Pr 

r 

wall BC 

plate BC 

N x 

x N z 

Nbl 

max(5a; ,5 Z ) 

Nu 

T f 

V 

lo 8 

1 

1.50 

SF 

NS 

512 x 256 

13 

0.68 

29.1 

600 

2.15 x 10 s 

1 

1.50 

SF 

NS 

1024 x 512 

23 

0.41 

36.1 

400 

4.46 x 10 8 

1 

1.50 

SF 

NS 

1024 x 512 

19 

0.55 

45.1 

400 

10 9 

1 

1.50 

SF 

NS 

1024 x 512 

16 

0.70 

56.0 

400 

2.15 x 10 9 

1 

1.50 

SF 

NS 

2048 x 1024 

30 

0.45 

66.1 

400 

4.64 x 10 9 

1 

1.50 

SF 

NS 

2048 x 1024 

25 

0.58 

83.7 

400 

10 10 

1 

1.50 

SF 

NS 

3072 x 1536 

36 

0.49 

101 

300 

10 8 

1 

2.00 

SF 

NS 

512 x 256 

9 

0.90 

29.7 

500 

2.15 x 10 8 

1 

2.00 

SF 

NS 

576 x 320 

11 

1.03 

36.4 

600 

4.46 x 10 8 

1 

2.00 

SF 

NS 

1024 x 512 

17 

0.74 

45.8 

600 

10 9 

1 

2.00 

SF 

NS 

1024 x 512 

15 

0.95 

56.8 

500 

2.15 x 10 9 

1 

2.00 

SF 

NS 

1152 x 576 

15 

1.07 

69.0 

500 

4.64 x 10 9 

1 

2.00 

SF 

NS 

1280 x 640 

15 

1.22 

81.9 

600 

10 10 

1 

2.00 

SF 

NS 

2048 

X 

1440 

41 

0.98 

101 

500 

10 8 

1 

1.00 

SF 

NS 

512 

X 

512 

18 

0.47 

33.4 

400 

2.15 x 10 s 

1 

1.00 

SF 

NS 

512 

X 

512 

17 

0.57 

35.4 

500 

4.46 x 10 8 

1 

1.00 

SF 

NS 

512 

X 

512 

19 

0.74 

44.5 

500 

10 9 

1 

1.00 

SF 

NS 

512 

X 

512 

11 

0.94 

54.6 

500 

2.15 x 10 9 

1 

1.00 

SF 

NS 

1024 

X 

1024 

29 

0.61 

70.1 

600 

4.46 x 10 9 

1 

1.00 

SF 

NS 

1024 

X 

1024 

23 

0.79 

92.9 

500 

10 10 

1 

1.00 

SF 

NS 

2048 

X 

1441 

35 

0.49 

103 

200 

2.15 x 10 10 

1 

1.00 

SF 

NS 

2048 

X 

1441 

28 

0.64 

135 

200 

4.64 x 10 10 

1 

1.00 

SF 

NS 

2048 

X 

1441 

22 

1.00 

180 

300 

10 11 

1 

1.00 

SF 

NS 

3240 

X 

3072 

59 

0.66 

217 

100 

10 8 

1 

0.33 

SF 

NS 

512 

X 

512 

16 

0.16 

37.8 

500 

2.15 x 10 s 

1 

0.33 

SF 

NS 

512 

X 

512 

13 

0.20 

46.5 

500 

4.46 x 10 8 

1 

0.33 

SF 

NS 

512 

X 

512 

11 

0.26 

57.0 

500 

10 9 

1 

0.33 

SF 

NS 

1024 

X 

1024 

35 

0.17 

73.5 

200 

2.15 x 10 9 

1 

0.33 

SF 

NS 

640 

X 

768 

18 

0.35 

97.8 

400 

4.46 x 10 9 

1 

0.33 

SF 

NS 

640 

X 

768 

15 

0.45 

125 

400 

10 10 

1 

0.33 

SF 

NS 

1024 

X 

1280 

27 

0.36 

154 

500 

2.15 x 10 10 

1 

0.33 

SF 

NS 

1441 

X 

1600 

31 

0.33 

199 

400 

4.64 x 10 10 

1 

0.33 

SF 

NS 

2048 

X 

3072 

66 

0.29 

240 

100 

10 11 

1 

0.33 

SF 

NS 

3072 

X 

3456 

69 

0.25 

280 

200 

10 8 

1 

0.30 

SF 

NS 

512 

X 

512 

17 

0.15 

38.5 

400 

10 8 

1 

0.35 

SF 

NS 

512 

X 

512 

17 

0.17 

37.3 

500 

10 8 

1 

0.40 

SF 

NS 

512 

X 

512 

18 

0.19 

36.1 

500 

10 8 

1 

0.45 

SF 

NS 

512 

X 

512 

18 

0.21 

35.0 

500 

10 8 

1 

0.50 

SF 

NS 

512 

X 

512 

19 

0.23 

34.0 

500 

10 8 

1 

0.55 

SF 

NS 

512 

X 

512 

19 

0.26 

33.1 

500 

10 8 

1 

0.60 

SF 

NS 

512 

X 

512 

20 

0.28 

32.0 

400 

10 8 

1 

0.65 

SF 

NS 

512 

X 

512 

20 

0.30 

31.5 

500 

10 8 

1 

0.70 

SF 

NS 

512 

X 

512 

21 

0.32 

31.0 

500 

10 8 

1 

0.75 

SF 

NS 

512 

X 

512 

21 

0.34 

30.5 

500 

10 8 

1 

0.80 

SF 

NS 

512 

X 

512 

21 

0.36 

29.6 

400 

10 8 

1 

0.85 

SF 

NS 

512 

X 

512 

21 

0.38 

29.5 

500 

10 8 

1 

0.90 

SF 

NS 

512 

X 

512 

22 

0.40 

28.5 

400 

10 8 

1 

0.95 

SF 

NS 

512 

X 

512 

22 

0.43 

28.8 

500 
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TABLE III. The columns from left to right indicate Ra, Pr, T, the wall BC, the plate BC, the resolution in horizontal and 
vertical direction N x x N z , the number of grid points in the thermal boundary layer #{n | n £ {1, 2, N z } A Z(n) < L/ (2Nu)} 
(where Z(n) is the vertical grid point distribution), the average Kolmogorov length scale in the flow compared to the largest 
grid length used somewhere in the grid max^, S z )/rj with r\/L — Pr 1//2 /[Ra(Nu — l)] 1 ^ 4 , the mean Nu over the full simulation 
and the averaging time in free fall time units r/. 


Ra 

Pr 

r 

wall BC 

plate BC 

N x x N z 

Nbl 

max(5 cc ,5 Z ) 

Nu 

T f 

V 

10 s 

1 

1.00 

NS 

NS 

256 x 256 

10 

1.12 

26.1 

300 

2.15 x 10 s 

1 

1.00 

NS 

NS 

512 x 512 

20 

0.76 

32.6 

300 

4.46 x 10 8 

1 

1.00 

NS 

NS 

512 x 512 

18 

1.00 

40.8 

400 

10 9 

1 

1.00 

NS 

NS 

512 x 512 

19 

1.33 

50.5 

800 

2.15 x 10 9 

1 

1.00 

NS 

NS 

1024 x 1440 

52 

0.85 

61.5 

400 

4.46 x 10 9 

1 

1.00 

NS 

NS 

1024 x 1440 

44 

1.08 

75.3 

400 

10 10 

1 

1.00 

NS 

NS 

1024 x 1440 

37 

1.45 

98.0 

200 

2.15 x 10 10 

1 

1.00 

NS 

NS 

2048 x 2048 

55 

0.94 

125 

200 

4.64 x 10 10 

1 

1.00 

NS 

NS 

3072 x 3072 

85 

0.80 

145 

100 

10 11 

1 

1.00 

NS 

NS 

3072 x 3072 

72 

1.04 

187 

200 

10 8 

1 

0.33 

NS 

NS 

256 x 256 

15 

0.36 

19.2 

1000 

2.15 x 10 s 

1 

0.33 

NS 

NS 

512 x 512 

13 

0.34 

23.2 

1000 

4.46 x 10 s 

1 

0.33 

NS 

NS 

512 x 512 

14 

0.38 

29.8 

900 

10 9 

1 

0.33 

NS 

NS 

512 x 512 

24 

0.76 

36.6 

300 

2.15 x 10 9 

1 

0.33 

NS 

NS 

512 x 576 

17 

0.49 

46.3 

200 

4.46 x 10 9 

1 

0.33 

NS 

NS 

512 x 1024 

33 

0.63 

59.3 

300 

10 10 

1 

0.33 

NS 

NS 

768 x 1440 

43 

0.59 

87.9 

200 

2.15 x 10 10 

1 

0.33 

NS 

NS 

768 x 1440 

34 

0.77 

120 

300 

4.64 x 10 10 

1 

0.33 

NS 

NS 

2048 x 2048 

26 

1.01 

172 

300 

10 11 

1 

0.33 

NS 

NS 

2048 x 3072 

63 

0.51 

200 

100 

10 7 

40 

2.00 

PD 

SF 

2048 x 1024 

41 

0.02 

37.4 

200 

10 8 

40 

2.00 

PD 

SF 

3072 x 2048 

56 

0.03 

87.9 

200 

10 9 

40 

2.00 

PD 

SF 

4096 x 2048 

40 

0.05 

168 

100 

10 6 

1 

0.50 

PD 

NS 

256 x 256 

65 

0.07 

2.61 

90000 

10 7 

1 

0.50 

PD 

NS 

512 x 512 

113 

0.07 

3.41 

300000 

10 8 

1 

0.50 

PD 

NS 

512 x 512 

96 

0.13 

4.32 

200000 

10 9 

1 

0.50 

PD 

NS 

1024 x 1024 

171 

0.13 

5.79 

20000 

10 11 

1 

0.50 

PD 

NS 

2048 x 3072 

105 

0.43 

95.4 

200 

10 12 

1 

0.50 

PD 

NS 

3072 x 4096 

73 

0.67 

282 

100 





